Aberration-like cusped focusing in the post-paraxial 
cn": Talbot effect 

o 

James D Ring 1 , Jari Lindberg 1 , Christopher J Howls 2 and 
Mark R Dennis 1 

. 1 H. H. Wills Physics Laboratory, University of Bristol, Tyndall Avenue, Bristol, BS8 

\Q i 1TL, UK. 

£N| ' 2 School of Mathematics, University of Southampton, Highfield, Southampton, S017 

1BJ, UK. 

^ . E-mail: james.ring@bristol.ac.uk 



O 



43 
Oh 



O 

(N 



Abstract. We present an analysis of self-imaging in a regime beyond the paraxial, 
where deviation from simple paraxial propagation causes apparent self-imaging 



C/3 . 

Q , aberrations. The resulting structures are examples of aberration without rays and 



^3 ' are described analytically using post-paraxial theory. They are shown to relate to, 

but surprisingly do not precisely replicate, a standard integral representation of a 
diffraction cusp. Beyond the Talbot effect, this result is significant as it illustrates that 
the effect of aberration - as manifested in the replacement of a perfect focus with a cusp- 
like pattern - can occur as a consequence of improving the paraxial approximation, 

T 1 

| rather than due to imperfections in the optical system. 

ON 
Q\ , 

I PACS numbers: 42.25.Fx, 42.15.Fr, 42.25.Hz 

t^j- \ The Talbot effect [U [2] is a well-known self-imaging phenomenon in classical 

optics, where a transversely periodic field reproduces itself at periodic distances in the 
propagation direction. As a natural consequence of paraxial wave propagation, it has 
also been observed in the diffraction of matter waves [3] and has an analogue in the 
revivals of quantum wavepackets [I]. Recently, the Talbot effect for surface plasmons 
[21 El [7J El [9] has attracted significant attention because of applications in areas such as 
nanolithography [TOj [H] . 

The perfect image revival of the classical Talbot effect can be understood in 
terms of Fresnel diffraction [2] and interpreted as an artefact of the corresponding 
paraxial approximation. Refining the propagation beyond the paraxial approach leads 
to aberration-like distortion of perfect self-imaging (T2J [13l [14] . Nonparaxial self-imaging 
has been studied previously for binary gratings using both scalar theory (TSJ [16] and full 
electromagnetic treatment [17] , and has particular importance in plasmonics where the 
transverse periodicity is comparable to the plasmon wavelength 0. Such nonparaxial 
self-images are not perfect revivals as they are in the paraxial case, and are only 
approximately periodic in the direction of propagation. 
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In this paper, we investigate a nonparaxial Talbot effect and explain that the 
resulting apparent aberration - the maximum intensity shifts and distorts to a cusp-like 
pattern - is not caused by imperfection in a lens or optical system (since the Talbot 
effect is lensless imaging), but rather as an actual improvement of the approximation, 
appropriate to our diffraction regime. Specifically, we invoke the fourth-order post- 
paraxial approximation [121 US] to investigate the Talbot-based focusing of scalar 
waves in a regime which is nonparaxial, but not extremely so: the transverse periodicity 
is of the order of some tens of wavelengths. This explains the appearance of diffraction 
cusps in the intensity patterns of surface plasmons calculated in Ref. [5]. We show that 
these aberration-like effects are not unique to plasmons and are connected to a standard 
integral representation of the diffraction pattern about a cusped caustic [18] despite the 
absence here of any actual ray caustic. Furthermore, we obtain a simple expression for 
the approximate location of the post-paraxial foci, which are shifted to be short of the 
paraxial Talbot distance. 

To demonstrate the physics of the post-paraxial regime, we employ a general two- 
dimensional (2D) formalism (transverse in x and propagating in z) which displays 
properties of both paraxial and nonparaxial regimes in different limits. Specifically, 
it consists of a transverse array (with period a in the x-direction) of monochromatic 2D 
point sources of wavelength A. The interference of the propagating waves in the x-z plane 
weaves a diffraction carpet, which is similar to that studied in the plasmon Talbot effect 
of [5]. The paraxial approximation applies when 7 = a/A is asymptotically large. In 
this case the Talbot distance 2a 2 /A is the spacing between planes of image revival. The 
initial pattern also revives at odd multiples of a 2 /A, albeit transversely shifted by a/2, 
leading to an alternate definition of the Talbot distance [15J of z^ = a 2 /X (used in this 
paper). In what follows, we use the dimensionless variables X = x/a and Z = z/zj?. In 
general, for any 7, the field at a distance Z, propagating according to the 2D Helmholtz 
equation, is given by (cf. Eq. (1) of [5]) 

N 

il>(X,Z)= exp 

n=—N 

where we assume negligible attenuation and consider only propagating waves, i.e., 
iV = [7] with square brackets denoting the integer part. Eq. [T]can be approximated by 
expanding 

y/1 - n 2 /-f 2 « 1 - n 2 /2 7 2 - n 4 /8 7 4 . (2) 
In the paraxial approximation to Eq. [H the quartic term in Eq. [2] is negligibly small, so 

TV 

i> p (X, Z)= ex P I 1 fanX " k^Z)] , (3) 

n=-N 

where we have neglected the overall phase factor of exp(27riz/ A). As shown in Fig. Ufa), 
the initial paraxial field revives for integer Z, since the second term in Eq. |3] will 
introduce a factor of ±1, and we recover the same field as for Z = 0. The factor of ±1 



2vri (nX + fZ^l-n 2 ^ 



(1) 
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Figure 1. Intensity distributions (7 — 18) for (a) the paraxial field of Eq. [3J (b) exact 
field of Eq. [T] (c) numerically simulated ID binary gold grating (of refractive index 
0.2 + 3.44i, thickness 100 nm and slit width to grating period ratio of 0.05) illuminated 
by a normally incident, TE-polarized, monochromatic (A = 632.8 nm) plane wave; and 
(d) the post-paraxial field as given by Eq. |4j A value of 7 = 18 is outside the region 
of validity for the paraxial approximation, since (a) is fundamentally different to (b). 
The fields (b) , (c) and (d) exhibit foci (indicated by the white box) which are aberrated 
compared to (a). 
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Figure 2. A comparison of foci symmetric about the centre of each image (7 = 18, 
M = 9) with colours corresponding to Fig. [TJ (a) paraxial focus of Eq. [3] (b) the exact 
field of Eq. [T] (c) numerically simulated focus according to the parameters of Fig.[TJc); 
(d) post-paraxial focus of Eq. [5] < 1); (e) the Pearcey function Eq. [TJ and (/) the 
incomplete Pearcey function (/i = term of Eq. [5]). 



determines whether or not the self-image is shifted by half a period. A post-paraxial 
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refinement is made by retaining the quartic term in Eq. [2] [121 EH EES] , giving 

N 

ip pp {X, Z)= ex P [i (2vmX - 7rn 2 Z - 7m 4 Z/4 7 2 )] . (4) 

n=-N 

Comparison between Eq. H] and Eq. [3] shows that the perfect Z-periodicity of the 
Talbot effect is an artefact of the paraxial approximation, and is destroyed in the post- 
paraxial regime by the inclusion of the quartic (and possibly higher terms) [15] . Figs. W[p) 
and W[d) show the exact and post-paraxial intensity distributions respectively, and it is 
immediately apparent that while the foci of these fields are similar to one another, they 
are very different, and less defined, from those of the paraxial case, (a). Comparing the 
foci of each regime in Fig. [1] it is clear that the perfect foci of the paraxial approximation 
of (a) have been perturbed in the post-paraxial regime to cusp-like distributions, evident 
in (b) and (d). Fig.[Tfc) shows a separate simulation for a binary grating, which has been 
calculated numerically using rigorous electromagnetic theory [T71 [T9] (parameters given 
in the caption), which very closely resembles the exact and post-paraxial fields in (b) 
and (d), including the cusp-like foci. While such a pattern is typical for aberration, its 
appearance here is because of a refinement of the usual Fresnel paraxial approximation. 
In fact, the inclusion of a quartic term destroying a perfect focus is exactly analogous 
to the spherical aberration of a parabolic mirror, though remarkably, in Fig. [1] it is 
improvement of the approximation and not any aberration in the system that causes 
the effect! 

The fields in Fig. [1] are those for an infinite array, and applying the paraxial 
approximation to such a system is not fully consistent [16], since at all Z there would 
be oblique contributions to the field from arbitrarily transversely distant sources - a 
fact that is at odds with the concept of paraxiality. Instead we model a finite array 
of N sources by truncating the sums of Eqs. [H, [3] and H] at ±M, where M = [N/2Z] 
[T5| [20] , embodying the walk-off effects of finite arrays [2] . Fig. [2] shows the various 
foci according to this truncation. Comparing the paraxial focus Fig. Wji) to the exact 
Fig. Wib)-, the aberration effect is clear. Fig. E^c) shows the cusp pattern from the 
diffraction grating of Fig. [H(c), but includes only diffraction orders up to ±M. There is 
a slight discrepancy between (6) and (c) because of the equal amplitudes of the planes 
waves in the superposition of Eq. [1] compared to varying amplitudes of the diffraction 
orders of the simulation of (c). The fields are nevertheless very similar and the cusp 
reveals the breakdown of the paraxial approximation. Below, we develop an analytic 
post-paraxial approximation to these fields, shown in Figs. W(d), W(d), from inclusion of 
the quartic term (and no higher), which is enough to capture the structure of the exact 
and simulated foci. Including higher terms improves the quantitative agreement, but 
without further qualitative changes. 

To derive a simple representation of the complicated interference pattern, we follow 
a similar technique to Ref. [15] but employ a finite Poisson sum formula [21] to account 
for the truncation at ±M of the spectrum. Applying this to ifr^X, Z) yields the main 
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result of this paper, 

^pp(^C) = g + — 2 — ^ / ds #( s ;^C,/4 (5) 

where f , C, A*) = ^(C)exp {-i [s 4 + crz(C)s 2 + cr x (£, /i, and 
/(n; £, C) = exp {ivr [2£n - C^ 2 - (T + Qn 4 /^ 2 ] } ; 

/»(C) = v^y/[*CF + C)]*; 

r(C) = M/h(C)] 

f 7x(e,/i,c) = 27r(e-/i)M0; 

ffz (C) = ^(OC- (6) 

In the above, we have set Z = T + £, where £ is a small shift from the T th Talbot 
distance, and £ = X — T/2. The terms f(±M;£,Q are end-point contributions from 
the truncation of the spectrum. We recognise that the integral in Eq. [5] is an incomplete 
version [22J of the Pearcey function (231 [23] (i- e - with a finite integration domain). The 
Pearcey function is given by 

/oo 
ds exp [i (s 4 + zs 2 + xs)] , (7) 
-oo 

which describes the generic diffraction pattern decorating a cusp caustic [HJ [23]. 

The sum in Eq. [5] represents an array of cusp-like foci in the x-direction, with 
transverse and longitudinal scalings given by ax and o~z- Although in optics it is 
usual to see the Pearcey function as decoration around a cusp caustic of rays that 
arise from a geometric optics approach [18], this is not necessary |25j. Since here we are 
outside the applicability of geometrical optics, the propagating function Eq. [1] cannot 
be approximated in terms of rays, and instead the Pearcey function has come from the 
quartic post-paraxial term in Eq. [31 Importantly, since the maximum of the Pearcey 
function is below the geometric focus, the maxima of the post-paraxial self-images are 
below the paraxial Talbot distance. This differs from the nonparaxial Talbot distance 
[H [T7] which is accurate when only a small number of diffraction orders contribute. 

Fig. W^d) shows the diffraction pattern near the Talbot distance (Z — 1) according 
to Eq. [5] for \fi\ < 1. The intensity of the Pearcey function, scaled to the size set by 
o~x and o~z, is shown in Fig. EJ^e). Fig. W^f) is the incomplete Pearcey function (the 
[i = term of Eq. [5]). Comparing (6) and (d), we see that although the sum in Eq. [5] 
has infinitely many terms, only the first few are needed to accurately describe the fine 
details of the field about the focus. The transverse fringes come from the truncation of 
the integral in Eq. and subsequently are absent from the Pearcey function, (e), and 
are just evident in (/). In all these images, it is clear that the maximum intensity is 
short of the paraxial Talbot distance. We now obtain an expression for the approximate 
location of the maximum intensity. 

As a first approximation we take only the \i = term of Eq. [5] and let r(Q go 
to infinity. This is equivalent to approximating the focal region near the first Talbot 
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Figure 3. The intensity of (a) the exact field, Eq. [T] (normalised by (2[7/2] + 1)) and 
(b) the post-paraxial field (Eq. [5j summing < 1) as a function of 7. The dashed 
line is the approximate location of maximum intensity according to Eq. |8l The dots 
are the actual maxima of the exact field for integer 7. The solid line is the Talbot 
distance. 



distance with the Pearcey function, Eq. [0 By using the numerically calculated, constant 
position of maximum intensity of the unsealed Pearcey function, Zp e w —2.199, we can 
employ the scalings of Eq. [6] to calculate the approximate location of the maximum 
intensity of the post-paraxial focus as a shift, Cpeak, from the T th Talbot distance. By 
setting <J Z (Cpeak) = Z Pc , we obtain 

Cpeak = (zl e + Zp^Zl + ien^T^j j 8tt 7 2 . (8) 

Fig. |3]^a) shows the intensity of the exact field along the symmetry axis of the cusp £ = 
for —0.1 < C < 0.05 as a function of increasing 7. Asymptotically large 7 gives rise to 
maximum intensity at the Talbot distance (C = in the paraxial limit), but for smaller 
7 it is located short of this. Fig. |^&) shows the post-paraxial approximation of Eq. [S] 
including only the terms < 1. The dashed curve is the location of the maximum 
intensity according to the approximation of Eq. |8j Note that this depends on T, such 
that consecutive foci will be shifted from their paraxial counterparts by an increasing 
amount. The dots indicate the maxima of the exact field (Eq. [1]) for integer 7, though 
they do not always correspond to well localised single peaks, as in the case of 7 = 12 
and 7 = 28. The intensity fringes that we observe in Fig. [3] are a result of interference 
between the \x = and the neighbouring fi = ±1 terms. Inclusion of more //-terms 
reconstructs the discontinuities of Fig. [3]^a) (which arise from the limits of the sum in 
Eq. [1] changing discretely with 7). 

To conclude, we have applied scalar post-paraxial theory to explain the observed 
cusp-like foci in the Talbot effect for the post-paraxial regime, which have replaced the 
point foci of paraxial theory. While the analytic form of these post-paraxial foci is related 
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to the standard description of diffraction about a cusp, there is a subtle distinction since 
it is the incomplete Pearcey function that occurs and not the standard Pearcey function. 
The finite integration domain of the incomplete Pearcey is of significant asymptotic 
importance when compared with the infinite domain of the complete Pearcey function 
of Eq. [7J This is clear when comparing intensity patterns - while both are cusp-like 
distributions, there is additional complex detail accounted for in the post-paraxial Talbot 
effect that is absent from the normal Pearcey function. Another important difference 
is that there is no ray caustic underpinning the structure of the post-paraxial cusp foci 
as there is with the complete Pearcey function. Despite their differences, the location 
of the maximum of intensity of the post-paraxial foci can be fairly well approximated 
using a scaled Pearcey function. 

It is curious to see that the effect of post-paraxiality on a focus appears to be 
that of aberration, when compared to the paraxial case. However, while the analytic 
form mimics that for aberration (hence the inclusion of the Pearcey function, albeit 
incomplete), it actually embodies an improvement of the paraxial approximation for 
the post-paraxial regime. The phenomenon is a universal effect that occurs as the 
parameters of a system (regardless of its exact realisation) move away from paraxiality. 

The results in this paper may plausibly find application in plasmonics, since it 
is a field in which the pertinent wavelengths and periodicities are comparable, and so 
nonparaxial effects are important. While the lengthscales considered here are slightly 
greater than those experimentally investigated so far, it would be both interesting and 
challenging to confirm the aberration- like effects of post-paraxiality on the Talbot effect. 
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